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ABSTRACT 

We discuss the possibility of generating the baryon asymmetry of the Universe when 
the temperature of the Universe is much below the electroweak scale. In our model 
the evaporation of primordial black holes or the decay of massive particles re-heats the 
surrounding plasma to temperatures above the electroweak transition temperature lead- 
ing to the restoration of electroweak symmetry locally. The symmetry is broken again 
spontaneously as the plasma cools and a baryon asymmetry is generated during the phase 
transition. This mechanism generates sufficient asymmetry for a second order electroweak 
phase transition. For a first order phase transition, sufficient asymmetry is generated if 
viscous effects slow down the heated plasma as it moves away from the black hole. In our 
scenario there is no wash-out of the asymmetry after the phase transition as the plasma 
rapidly cools to lower temperatures thereby shutting off the sphaleron processes. 



I. INTRODUCTION 



In recent years much effort has been devoted to formulating a mechanism for baryo- 
genesis at the electroweak phase transition [1J . Electroweak baryogenesis is indeed a very 
exciting possibility as here one is working within the framework of a theory (or its exten- 
sions) which is reasonably well understood and the energy scale involved is accessible in 
laboratory experiments. However the requirement of a strong first order phase transition, 
large CP violation and a small Higgs mass to ensure that the created baryon asymme- 
try is not washed away after the phase transition places stringent constraints on these 
models and practically rules out electroweak baryogenesis in the context of the Standard 
Model. In this paper we discuss an alternative scenario for electroweak baryogenesis. Here 
the baryon asymmetry is created at temperatures much below the electroweak transition 
temperature during the evaporation of primordial black holes. When a black hole is evap- 
orating it heats up the plasma around it to a temperature much higher than the ambient 
temperature for a short time. This can also happen due to the decay of massive particles. 
For appropriate black hole masses (or, particle masses) the temperature of the hot region 
rises above the electroweak transition temperature T ew and the electroweak symmetry is 
restored locally. Due to the transfer of energy out of this region the hot region will cool 
and the temperature will fall below T ew . Thus in these hot regions the electroweak phase 
transition occurs again. Baryon asymmetry is then generated in these hot regions during 
this phase transition. 

One motivation for our model is that a first order transition at a very early stage in 
the Universe, or density fluctuations in general, seem to inevitably lead to the formation 
of primordial black holes. Evaporation of these black holes will lead to local heating of the 
plasma and hence to the possibility of electroweak baryogenesis for appropriate black hole 
masses. Also, scalars found in string theory can decay after the electroweak transition 
and dilute any baryon asymmetry created earlier ||. But in our scenario, we envision 
that the particles that dilute the baryon asymmetry during their decay can also recreate 
the baryon asymmetry after their decay. 

The basic mechanism of our model can be realized either by particles emitted in the 
evaporation of a primordial black hole or by particles constituting the decay products of a 
heavy particle. In both cases, it is the energy of these emitted particles which thermalizes 
and leads to local heating of the surrounding plasma. For the black hole case the heated 
region will be spherically symmetric while for the particle decay case the region will be 
collimated. In this paper we will only discuss the black hole case. The particle decay case 
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is much more complicated. We will briefly comment on this case at the end of the paper 
but it is not clear to us at this stage if in that case one can get sufficient heating of the 
local regions. We hope to discuss the particle decay case in a future work. 

A desirable feature of our scenario is that it works for a second order transition. In 
the standard cosmological scenario where the Universe cools due to Hubble expansion, 
a second order phase transition can not produce sufficient baryon asymmetry as the 
baryon asymmetry produced is proportional to 0/0 ~ H, where is the Higgs field, and 
this quantity is too smallQ In our scenario the heated plasma cools much faster as it 
moves away from the black hole and 0/0 is large enough to produce the observed baryon 
asymmetry. For a first order phase transition, the asymmetry generated is suppressed 
by the Lorentz factor associated with the relativistic speeds at which the heated plasma 
moves away from the black hole. However if viscous effects in the plasma slow down the 
outward moving plasma this suppression will be decreased. In our scenario, since the 
heated plasma cools rapidly, the temperature after baryogenesis quickly becomes much 
smaller than the electroweak scale preventing any significant wash-out of the created 
asymmetry after the phase transition. 

We divide the calculation of the baryon asymmetry in our model into two main steps. 
In Section 2, we estimate the volume of the plasma in which the electroweak symmetry 
is restored and which is relevant for baryogenesis in our scenario. In Section 3, we esti- 
mate the baryon asymmetry generated as these regions cool and undergo the electroweak 
transition. Then in Section 4 we point out why the upper bound on the Higgs mass does 
not exist in our model. We conclude in Section 5 with a brief discussion of how inhomo- 
geneities in the baryon number distribution can be smoothed out before nucleosynthesis. 

We briefly comment on the formation of primordial black holes in the early Universe. 
In pre-inflationary days it was argued that the initial spectrum of density inhomogeneities, 
invoked to explain the observed structure in the Universe, could also give rise to primor- 
dial black holes It was also argued that white holes would be unstable and would 
convert to black holes 0. In first order phase transitions, particularly in first order in- 
flationary scenarios, primordial black holes can be produced by collapsing regions of false 
vacuum trapped between bubbles of the true vacuum Primordial black holes can 



* However models of electroweak baryogenesis with topological defects are insensitive to the 
order of the electroweak phase transition and can create the observed asymmetry for a second 
order phase transition Q. 
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also be formed via the gravitational instability of inhomogeneities formed during bubble 
wall collisions in first order inflation H-ff^], by large amplitude density perturbations pro- 
duced due to fluctuations in the inflaton field |TB| , by shrinking cosmic string loops |H| . [To'] 



and by expanding topological defects produced during inflation [Tj| . As our primary mo- 
tivation is to demonstrate a qualitatively different scenario for implementing electroweak 
baryogenesis, we will not go in to the specific mechanism which could give rise to the 
formation of black holes of required masses. However we point out that the formation of 
black holes which evaporate below the electroweak scale has been discussed in [fT7| . 



In ref. [18] a scenario which is similar to ours has been proposed. We comment later 



on this work. The presence of a hot plasma surrounding primordial black holes has also 



been discussed in refs. ||i9|-f22f . However, we point out that in these references the plasma 
consists of the black hole radiation and particles produced by the interaction of the black 
hole radiation with itself. In contrast, in our scenario we consider the interaction of the 
Hawking radiation with the ambient plasma surrounding the black hole, and the subse- 
quent heating up of this plasma by the Hawking radiation. The possibility of obtaining 
the baryon asymmetry of the universe in the plasma around the black hole is mentioned 
in refs. [0,0] . 

II. FORMATION OF HOT REGIONS 

A black hole evaporates by emitting Hawking radiation with an associated temperature 

t - M2pi m 

where M b h is the mass of black hole and Mp\ = 1.2 x 10 19 GeV is the Planck mass. We 
use natural units with H — c — 1. The rate of loss of mass by the evaporating black hole 
is given by 



dM bh aM^ 
dt ~ ~ Ml 



(2) 



L bh 

Here, a accounts for the scattering of emitted particles by the curvature and depends 
on T bh . For different values of T bh values of a have been tabulated in p3[. For T bh = 



1, 200 and 10 15 MeV, the corresponding values of a are 3.6 x 10 -4 , 2.3 x 10" 3 and 4.5 
xlO -3 respectively. As we will see later, the relevant value of T b h for us will be higher 
than about 10 6 GeV. Therefore we shall set a to be 3 x 10 -3 . 

The lifetime Tbh of the black hole can be obtained by integrating eqn.(2). We get 
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r bh ~ 10 2 Mp; 4 M 3 



(3) 



where M is the initial mass of the black hole. Eqn.(2) implies that very little energy 
is emitted until time of the order of which is when most of the energy of the black 
hole gets emitted. Thus for a black hole formed early in the Universe, it is reasonable to 
assume that the black hole essentially evaporates only when the age of the Universe is of 
order r bh . 

Let us assume that black holes formed in the early Universe have masses so that they 
evaporate when the temperature of the Universe is Ty < T ew . The age of the Universe ty 
when its temperature is Ty is 

t v ~ 0.3g^ 2 M Pl Ty 2 (4) 

Here g* is the number of degrees of freedom relevant at temperature Ty. We are 
interested in black holes that decay after the electroweak phase transition but before the 
onset of nucleosynthesis at Ty — 1 MeV. For concreteness, we shall consider black holes 
decaying at Ty equal to 1 GeVand 10 GeV. For Ty in the range 1-100 GeV, is equal 
p4| to 100. From eqn.(3) we can then calculate M so that r b h = ty. We get (by using g* 
= 100), 

M = 0.07Mpf Ty 2/3 (5) 

The temperature of this black hole is 

T bh = O.QM l P { 3 Ty /3 . (6) 

For Ty = 1 GeV we get M = 4 x 10 n M Pi , T bh = 1 x 10 6 GeV and r bh = 5 x 
10 17 GeV -1 = 3 x 10~ 7 s. The picture then is that these black holes emit particles with 
energies roughly equal to 10 6 GeV into the background plasma which is at a temperature 
Ty (~ 1 GeV) to start with. These 10 6 GeV particles will scatter with the particles 
in the background plasma and will heat it up through their energy loss. For the black 
hole masses considered here only elementary particles will be emitted, such as quarks, 
gluons, photons, leptons, etc. (Emission of quarks and gluons by black holes is a non- 
trivial process as discussed in P5| . We will not worry about those details here.) For 
Ty = 10 GeV, M = 8 x 10 10 M„, T bh = 6 x 10 6 GeV and r tfl =4x 10 15 GeV" 1 = 3 x 10~ 9 
s. 

Obtaining the temperature profile outside a black hole radiating into an ambient 
plasma is non-trivial. We first show below how energy from the black hole is deposited 
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in the surrounding plasma thereby heating it up. It is necessary to include heat trans- 
fer to allow for the energy deposited close to the black hole to move out. We first use 
diffusion equations from stellar physics but then raise a concern that one does not have 
hydrostatic equilibrium in the plasma surrounding the black hole, unlike in stellar atmo- 
spheres. Therefore the plasma also moves out due to the pressure gradient. We model 
this outflow of the plasma using shells. However a more complicated numerical study may 
be necessary to understand the simultaneous heating of the plasma and energy transfer 
due to diffusion and bulk flow. 

The energy loss of particles traversing a region of quark-gluon plasma has been dis- 
cussed extensively in the literature. The energy loss per unit distance for an energetic 
quark with energy E traversing a quark-gluon plasma at temperature T is given byf\ f26fl 



where we have taken a s ~ 0.1. For the black hole of mass 4 x 10 n Mp/, setting T equal 
to the ambient temperature Tjj = 1 GeV and AE = 10 6 GeV, we get Ax ~ 10 6 GeV _1 , 
if we approximate E(x) as T^. One could argue that Aa; is the distance in which the 
energy released by the black hole is thermalised as by Ax an emitted particle has attained 
an energy comparable to the thermal energy of the particles in the surrounding plasma 
(loosely speaking, it has thermalised). Furthermore, it takes many collisions before the 
radiated particle comes to the same energy as the background, in which time the heat 
transferred along the way would also thermalise due to collisions. If one now equates 
Mbh = |^-^*T 4 ^ L r^ where the stopping distance r st = Ax then T = 400 GeV. However we 
are here grossly assuming that all the emitted energetic particles go through a plasma at 
1 GeVand that all the black hole energy is thermalised to a uniform temperature. Instead 
let us consider the black hole evaporation in steps and equate COIM^ with the increase in 
energy of the ambient plasma within r st . The temperature of the plasma rises to 100 GeV. 
The next 1% of the black hole mass released as Hawking radiation will now see an ambient 
plasma at 100 GeVand hence the new stopping distance r' st will be 200 GeV -1 . Thus a 
smaller region than before is heated to an even higher temperature. This continual process 
can give rise to a temperature gradient in the plasma surrounding the plasma. The above 
also shows that naively equating the energy emitted with p^-r^ t , as we did above, will not 



tThis expression was derived for energies lower than 10 6 GeV but since a s is less for higher E 
the perturbative derivation of the energy loss is even more accurate for our case. 




(7) 
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be correct. Though the plasma within r st rises to T ew even with the energy of 1% of the 
black hole we do not use this to estimate the volume in which the symmetry is restored 
as it is not clear what fraction of the black hole mass, emitted as radiation, sees a plasma 
at the initial ambient temperature of 1 GeV. 

To obtain the temperature profile outside the black hole one must include heat transfer 
from the regions close to the black hole to outer regions. (In fact, if we ignore heat 
transfer the temperature of the plasma close to the black hole keeps rising until the plasma 
temperature becomes equal to the black hole temperature. After this the black hole can 
not evaporate any further until the energy deposited close to the black hole moves out.) 
To include the effects of diffusion we shall follow the calculation of the temperature profile 
of stellar interiors. In stars one has a hot stellar core in which nuclear reactions occur, and 
a stellar atmosphere surrounding the core. Energy moves through the stellar atmosphere 
via diffusion due to a thermal gradient. We have an analogous situation with a black hole 
radiating into a plasma around it and we shall assume that energy radiated by the black 
hole is quickly thermalised, so at any given time one has a hot black hole radiating into 
a plasma in which there is a thermal gradient. Energy diffuses from the hotter regions 
closer to the black hole to regions further away due to the thermal gradient. Under these 
assumptions we will use eqn. 5.11 of |27J (hereafter referred to as KW), namely, 

lbnac r A l 4 

where c = 1 is the speed of light, a = (ii 2 /30)g il: . is the radiative cross-section per 
unit mass averaged over frequency times the energy density and equals an =l/mean free 
path, where a is the cross section for scattering and n is the number density of particles. 
The local luminosity l(r) is defined as the net energy passing outward through a sphere 
of radius r per second (see Section 4.2 of KW). 

To integrate dT/dr, we need an and l{r). We use 

an = (1/9)(2) 3(4^) X ( 3 /4)(l-2/vr 2 )^T 3 • (9) 

The above expression for a is the cross section for quark- ant iquark scattering (see eqn. 
6.33 of ||28|| ). The factors of (l/9)(2) come from averaging over color and tracing over 
QCD generators, respectively. We have set the Mandelstam variable s in the formula to 
equal (2T) 2 , i.e., for incoming particles of equal (and opposite) momenta T. Then we find 
that an = 2 x 10~ A g*T = 0.02T, where we have used a s = 0.1. 

In stars l(r) can be a complicated function, depending on the distribution of sources 
and sinks of energy. Since we are interested in the steady state situation we shall take 
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the luminosity to be independent of r and equal to the energy radiated per unit time by 
the black hole. Therefore 



l{r) = -dM bh /dt = aM Pl /M^ h . (10) 

Now one can integrate dT/dr, with the boundary condition that T = T' at r = r' . We 
shall initially take T' to be the ambient temperature Tjj. A priori one does not know r'. 
Due to scattering of radiation in the plasma r' should be less than cr^. We shall fix r' by 
requiring that the increase in the energy of the plasma for r < r' is equal to the radiated 
mass of the black hole. We get 



(9/l 6vrac )__^(0.02) 



1 1 



(11) 



r r 

We shall apply this solution for r > QR, where R = 2Mj ) h/Mp l is the Schwarzschild radius 
of the black hole. 6R is the radius of the thermal atmosphere around the black hole |2"9"| 
beyond which curvature effects of the black hole are negligible. 

Our assumption of l(r) independent of r implies a steady state profile. However our 
choice of the boundary condition that the increase in energy of the plasma within r' equals 
the energy radiated by the black hole implies that all the black hole energy is used to 
create the temperature profile. More likely, the energy radiated from the black hole will 
first create a profile and then the remaining energy of the black hole will be transmitted 
through the plasma without further heating the plasma. This implies that the increase 
in energy of the plasma will be less than the radiated mass of the black hole. Therefore 
T(r) for any r will be less than what we obtain and our estimates of rioo in the diffusion 
picture will be an upper bound. 

We start with a black hole of mass M bh = 4x 10 n M Pl and lifetime 5 x 10 17 GeV -1 which 
decays when the temperature of the universe is 1 GeV. Using the temperature profile in 
eqn.(|Tl]) we apply the energy constraint mentioned above to (numerically) obtain r' = 
2 x 10 9 GeV -1 , and thus an expression for T{r). While applying the energy constraint we 
assume that only 90% of the black hole mass is radiated away. Using T(r) we then obtain 
riooj the radius below which the temperature is greater than 100 GeV, to be 4 x 10 2 GeV -1 . 

Now we envision that the black hole has reduced in mass by 90% and so we have a 
black hole of mass 4 x 10 10 Mp; and lifetime r = 5 x 10 14 GeV -1 . We once again use the 
above equation and find that the energy constraint for a black hole that radiates 90% of 
its mass is consistent with r'= lx 10 8 GeV -1 and T(r') = 2 GeV. (We vary r' and T(r') so 
that we satisfy the energy condition for a value of r' such that T(r') is consistent with the 
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temperature profile for M bh = 4 x 10 11 M Pi . Also, the energy in the region r < r' obtained 
from the temperature profile of the black hole at M bh = 4 x 10 11 Mp; is much smaller 
than the energy deposited in this region by the black hole after its mass has reduced to 
4 x 10 10 M«.) We now find that r 100 = 3 x 10 4 GeV -1 . 

Now let the mass of the black hole be 4 x 10 9 M Pl . Its lifetime is 5 x 10 11 GeV -1 . We 
repeat the procedure above and obtain a consistent temperature profile with T' = 30 GeV 
at r' = 1 x 10 6 GeV -1 . Again the energy deposited greatly exceeds the energy already in 
the plasma for r < r'. We find that r 100 = 8 x 10 5 GeV -1 . Since r 100 ~ r' we can not get 
a larger r wo as the black hole mass continues to reduce further. Thus for a black hole of 
mass 4 x 10 11 M Pi radiating into a plasma (initially) at a temperature of 1 GeV we obtain 
r 100 = 8 x 10 5 GeV -1 . 

Had we started with an initial ambient temperature of 10 GeV, and a black hole of 
mass 8 x 10 10 M Pi , temperature 6 x 10 6 GeV and lifetime 4 x 10 15 GeV -1 or 3 x 10~ 9 s, we 
would first obtain r' = 8x 10 7 GeV -1 for V = 10 GeV and r 100 = 8 x 10 3 GeV" 1 . Allowing 
the black hole to further radiate with a mass of 8 x 10 9 GeV we find r' = 5x 10 6 GeV with 
V = 10 GeV. Then we obtain r 100 of 6 x 10 5 GeV -1 . 

If one calculates the energy radiated per unit time into the black hole by the plasma 
at r = 6R, P = a(T{6R))Mn(6R) 2 , where a = (1/4)(tt 2 /30)^ is the Stefan-Boltzmann 
constant modified for more than the two photonic degrees of freedom, one finds that this 
is always a factor of 10 3 less than —dM bh /dt\ vacuum = aM pl /M 2 h . Therefore we can ignore 
any decrease in the net luminosity due to absorption by the black hole. 

Though we have used the stellar diffusion equations to obtain the temperature profile 
there is one important difference between our situation and that of stellar atmospheres, 
namely, the existence of hydrostatic equilibrium. While gravity in stars provides hydro- 
static equilibrium we show below that pressure gradients in the plasma will overcome the 
gravitational attraction of the black hole and the plasma will move out. This is to be 
expected as these black holes have much less mass and plasma temperatures are much 
larger as compared to stellar situations. 

The gravitational force (per unit volume) on an element of the plasma of density p and 
at a distance r from the black hole is (p + P)g, where g = GnM/t 2 is the gravitational 
acceleration and P = |p is the pressure at a distance r from the black hole. M here 
includes the black hole mass and the energy plus pressure of the plasma within the sphere 
at radius r. One can show that the contribution of the plasma to M is negligible. Therefore 
the gravitational force (per unit volume) is 
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_4tt 2 _ T*M bh 
h ~ 90 9 * M 2 Pl r 2 ■ { ' 

The outward force (per unit volume) on the plasma element is given by VP = dP/dr. 
Now P = ^g^T 4 . Therefore 

An 2 

dP/dr = —g,T 3 dT/dr . (13) 
90 

For Mf,h = 4 x 10 u Mpi we find that outward force due to the pressure gradient is much 
larger than the gravitational attractive force for regions where T < 1 x 10 5 GeV. At r 1Q0 
dPj dr is 9 orders of magnitude larger than the gravitational force. For M^h = 8 x 10 10 Mpi 
gravity becomes sub-dominant at T = 7 x 10 5 GeV and the disparity between the pressure 
gradient and the gravitational force at r wo is even larger. Thus one should include effects 
of bulk flow in the modeling of the plasma outside the black hole. 

We now calculate the distance at which heat transfer becomes dominated by bulk flow 
rather than diffusion. As we have shown, beyond some distance r a hydrostatic equilibrium 
is not maintained and the plasma acquires some outward bulk velocity. However bulk 
plasma motion will not be the dominant mode of heat transfer till the plasma velocity 
increases sufficiently. Let r c be the distance at which heat transfer due to bulk flow 
becomes dominant. 



Starting from the Euler equation for a relativistic fluid [ 30 1 we obtain the following 



equation for steady-state situations (i.e. after setting time derivatives to zero) 

vrfdvjdr = -l/(p + P)dP/dr , (14) 



where 7 = 1/\A — u 2 and c = 1. Using the expression for the pressure of a relativistic 
fluid, P — (1/3) we get 

v/(l - v 2 )dv/dr = -(\/T)(dT/dr) . (15) 

We would like to integrate the above equation with the range of integration for v 
and r being to v c and r a to r c respectively. Since heat transfer in this region is 
still dominated by diffusion we use dT/dr = —B/(r 2 T 2 ) from eq. (|8|) above, where 
B = (90/(1600tt 3 ))0.06 x 10 _3 (M^/M 6 2 J. We obtain the expression for the tempera- 
ture profile T(r) as 

r 3 = T 3 + 3£(l/r-l/r ), (16) 

where To = rjTbh and ro = QR. Since we are trying to ascertain the distance r c at which the 
diffusion approximation breaks down, we can not impose a boundary condition at large r, 
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i.e., at r' as we have done earlier. However, our ignorance of the plasma temperature at 
small r = 6R leads us to introduce the parameter 77. 77 is then a measure of how different 
the temperature at r = 6R is from the black hole temperature T^. To know the precise 
value of 77 would require a complicated numerical simulation of the processes occurring 
close to the black hole. Therefore we instead solve for r c for different values of 77. For a 
black hole of mass 4 x 10 ll Mpi, if we substitute values of 77 > 0.12 in eq. (|16D then even at 
r = 00 the temperature is greater than T a = 1 x 10 5 GeV. Since hydrostatic equilibrium 
breaks down for temperatures below this value this implies that hydrostatic equilibrium is 
maintained for all r and that the diffusion approximation is valid everywhere. However the 
large value of T at infinity for these values of 77 implies that the entire universe is heated 
up by a single black hole which violates energy conservation and shows that 77 > 0.12 
should not be used. For 77 = 0.10, r a is 6R = 4 x 10 -7 GeV -1 and so for smaller values 
of 77 one also takes r a equal to 6R. For values of 77 between 0.10 and 0.12 r a is greater 
than QR and is obtained from eq. QTJ) to be 3.0 x 10 -10 (1.8 x 10 -3 - rj 3 )- 1 GeV -1 . For a 
black hole of mass 8 x 10 10 GeV values of 77 > 0.13 violate energy conservation while for 
77 < 0.12 we take r a as 6R = 8 x 10 -8 GeV -1 . For values of 77 between 0.12 and 0.13 r a is 
3.4 x 10 -11 (2.0 x 10 -3 - 77 3 ) -1 GeV -1 . 
v c is then obtained from eq. flTo] ) as 

v\ = 1 - exp(- £ 2B/(r 2 T 3 )dr^j . (17) 

r c can be obtained by equating the flux due to bulk flow at r c , J(r c ) = 7 2 (p + P)(r c )v c , 
with the total flux at r c , i.e., 

7 2 ^?*T C V = Ij^rl) (18) 

where T c is the temperature at r c and / is the luminosity at r c , which we take to be 
—dMbh/dt. This finally gives 

7 2 ^T C V> C = {aMU)l{^Ml) (19) 

By substituting the expressions for v c and T c obtained in terms of r c , it is in principle 
possible to obtain r c . However the above equation cannot be solved analytically. Instead 
we numerically estimate r c by plotting (LHS-RHS) of eqn.(|T9D as a function of r c to 
ascertain at what point the graph crosses the r c -axis. Substituting this value into eqs. ( |i~7D 
and (|16D gives v c and T c respectively. 
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As mentioned earlier, we present results for certain values of 77. For a black hole of mass 
Mbh = 4 x 10 11 Mp; evaporating when the ambient temperature of the Universe is Tjj = 

1 GeV we find r c = 9 x 1(T 6 GeV -1 and 3 x KT 6 GeV" 1 for 77 = 0.10 and 0.08 respectively. 
The corresponding values of v c are 0.5 and 0.8 respectively. The corresponding values of T c 
are 1 x 10 5 GeV and 6 x 10 4 GeV respectively. For a black hole of mass M bh = 8 x 10 10 M P i 
evaporating when the ambient temperature of the Universe is Ty = 10 GeV we find r c = 

2 x 10~ 6 GeV -1 and 5 x 10~ 6 for 77 = 0.10 and 0.08 respectively. The corresponding values 
of v c are 0.5 and 0.8 respectively. The corresponding values of T c are 5 x 10 5 GeV and 

3 x 10 5 GeV respectively. In all the above cases the temperature of the plasma at r c is 
much greater than 100 GeVindicating that r c is smaller than r 1Q0 . This shows that bulk 
flow dominates over diffusive transfer of energy in the plasma before the temperature 
of the plasma drops down to the electroweak transition temperature, thereby calling in 
question the validity of the diffusion picture. 

To study the outward motion of the plasma for r > r c we shall consider the plasma as 
consisting of outgoing shells of infinitesimal thickness dr (in the frame of the black hole). 
As these shells move out, they cool down due to the expansion of the plasma. Eventually 
these shells will reach a distance r wo where the temperature (as defined in the frame of the 
plasma) becomes smaller than 100 GeV. At this time the electroweak transition will occur 
in these shells and, as we will discuss later, baryogenesis can take place there. Rewriting 
eqn. (|l5|) and integrating we get 



which gives us a relation between the velocity of the outward moving plasma shell and 
the temperature of the plasma. 

Since the temperature decreases with increasing r, the velocity of the plasma shell will 
keep increasing and will eventually become very close to the speed of light. The Lorentz 
factor is given by 



For Tu — 1 GeV the Lorentz factor 7 100 at r 100 is 1000 for 77 = 0.10 and 0.08. For 
Tjj = 10 GeV, 7100 is 6000 and 5000 for 77 = 0.10 and 0.08 respectively. Below we shall 
take Vioo, the velocity of the plasma shell at r 10 o, to be 1. Later we shall briefly mention 
the possible role of viscosity in slowing down the plasma and its consequences. 

Assuming a steady state situation l(r) is independent of r and equal to —dM^h/dt. 
Therefore for r > r c 



v 2 = l-(T/T c ) 2 (l-v 2 ), 



(20) 
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(21) 
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Above T is the temperature of the plasma shell at a distance r. For a black hole of mass 
4 x 10 n M Pl , r 100 is 77 -1 GeV -1 and it is 3O7- 1 GeV -1 for a black hole of mass 8 x 10 10 M Pi . 

The volume relevant for baryogenesis must be calculated in the rest frame of the 
plasma since the sphaleron rates have been calculated in this frame. Consider a plasma 
shell of radius r wo in the black hole frame. There is no frame in which the entire shell is 
at rest. Instead one can divide the shell into plasma elements of area dflrf 0Q and thickness 
dr. In the rest frame of such a plasma element the volume element is dQrl 00 dr', where 
dr' = ^ydr is the shell thickness in the frame of the plasma element. Therefore the volume 
relevant for baryogenesis contributed by one shell is 

dV = 4nrl 00 dr' = ^r 2 w ^ im dr = 4:irr 2 00 j 100 v wo dt . (23) 

Plasma shells keep moving across r wo for a duration of the order of the lifetime r of the 
black hole. Thus the total volume relevant for baryogenesis is then given by 

V = 47rri 00 7i 00 vioo7- = ^1wo 1t / ( 4 Pioo) , (24) 

where we have used eq. ( p2|) for the second equality. Since I = —dM bh /dt = aM pl /M 2 h 
we may also express the volume as 

V = 7 x KT 11 GeV- 4 7i"oo M ^ = 5x 1(T 12 GeV^^^Mpf /T^ /3 (25) 

where we have used eq. (0) for the last equality. For a black hole of mass 4 x 10 11 Mp; 
radiating into a plasma initially at 1 GeV, V is 3 x 10 20 7{oq GeV~ 3 ; for a black hole of 
mass 8 x 10 10 Mpi and an initial ambient temperature of 10 GeV, V is 7 x 10 19 7 10 o 

Electroweak baryogenesis takes place via sphaleron processes with sphaleron size ~ 
{awT ew )~ l ~ (3 GeV) -1 . Since the thickness of the shell in the plasma frame for a 
temperature difference of 1 GeV,i.e., r' gg — r' 1Q0 , is proportional to a Lorentz factor, which 
compensates for the 7 factor in rioo, one can easily fit the sphaleron in the plasma shell 
at noo- 

III. ESTIMATION OF THE BARYON ASYMMETRY 

Let us assume that the number of excess baryons created around each black hole is 
Nb- If the number density of black holes is then the number density of baryon 
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excess produced by these black holes is ub = NsUbh- To estimate n^u we assume that the 
net energy density of black holes does not dominate the energy density of the Universe at 
the time when black holes decay. (With this assumption we are assured that black holes 
never dominate the dynamics of expansion of the Universe so the overall picture of the 
evolution of the Universe remains unchanged.) This gives the following expression for nbh 
(using eqn.(5)). 

n bh ^5x 10 2 Ti 4/3 M p f /3 (26) 

For Tu = 1 GeV, %<8x 1(T 30 GeV 3 , while for T v = 10 GeV, %<4x 1(T 25 GeV 3 . We 
shall use the maximum allowed value of n^h below. 

The entropy density of the Universe (which will remain essentially unchanged even 
after all the black holes evaporate, since the black holes do not dominate the energy 
density of the Universe) is given by 

s = %9*n ~ 40T 3 (27) 

We now proceed to estimate the baryon asymmetry produced in the hot regions by 
models of baryogenesis conventionally used at the electroweak scale. In the following we 
shall use V to refer to the hot region where the temperature rises to at least T ew as well 
as to the volume of this region. Let us first assume that the electroweak transition is of 
first order. Before the evaporation of the black hole the Universe is in a state with the 
expectation value of (ft at a value v. When the regions of volume V heat up, the expectation 
value of in V will evolve to the value 0. (Later we consider the possibility that <fi may 
not roll to the minimum due to insufficient heating.) As the plasma in V cools further, a 
local minimum at a non-zero value of appears in the effective potential. For T < T ew , 
this minimum becomes the true minimum of the potential. Bubbles of true vacuum now 
appear in the hot region V as the field tunnels from to v. As these bubbles grow, they 
either produce a baryon asymmetry in the outer regions of a thick wall by "spontaneous 
baryogenesis" |TiJ or by the asymmetric reflection of particles in the symmetric phase from 
the surface of a thin bubble wall |j32| . In the former case, the non-zero value of the time 
derivative of a field coupled to the baryon number acts as a chemical potential for baryon 
number, thereby biasing B+L violating sphaleron processes in the wall to produce more 
baryons than anti-baryons. In the latter case, the asymmetry created in some quantum 
number not orthogonal to baryon number due to reflections off the bubble wall is converted 
into a baryon asymmetry by sphaleron processes. This baryon asymmetry is frozen after 
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it passes through the wall. A more sophisticated treatment of CP violation, particle 
transport and diffusion in the context of electroweak baryogenesis is discussed in refs. 
|53|-57[]. A baryon asymmetry can also be produced at the electroweak phase transition 
through the CP asymmetric interaction of the bubble wall with fluctuations in the winding 
number of gauge-Higgs fields configurations [^,^] and through the interference between 
electroweak sphaleron-induced baryon number violating processes and QCD sphaleron- 
induced CP violating processes ]4(J. In the standard electroweak baryogenesis scenario 
baryogenesis occurs in the entire Universe as the electroweak bubbles sweep through 
the Universe during the phase transition. In our case, our hot regions cover a volume 
fraction of Vn^ = 3 x 10~ 9 7 _1 GeV _4 T^. However, the entropy density of the Universe 
at a temperature of Ty is also lower by a factor of 10 6 GeV 3 /T^ as compared to the 
entropy density at the electroweak phase transition. The ratio of the asymmetry that 
we obtain to that in the standard scenario is then volume fraction x entropy factor 
= 3 x KTV^c/GeV- 1 . 

The Standard Model does not have sufficient CP violation to create the observed 



baryon asymmetry of the Universe. Therefore we consider the results of ref. [36] for 
baryogenesis in supersymmetric models. With CP violating phases of the order of 10~ 4 
they produce enough asymmetry for a first order electroweak phase transition. However, 
even with a larger value of the CP violating parameters the asymmetry that we obtain 
above for a first order electroweak phase transition will be insufficient because of the large 
value of 7. 

However above we have completely ignored the effects of viscosity. Viscous effects can 
be of two types- bulk and shear []41|j . The bulk viscosity coefficient for a (quark-gluon) 
relativistic plasma is approximately |42|]. However there is a contribution to the shear 
viscosity that is proportional to dv/dr, which will slow down the plasma as it moves out. 
The relativistic hydrodynamic equations with viscosity are very complicated and can not 
be solved easily. If, however, the viscous effects decrease the plasma velocity at r wo to 
even 0.9c, this will imply that 7100 is 2. Then, with a value of CP violation ~ 0.1 or 
0.01, one can obtain sufficient asymmetry in our scenario for black holes radiating into an 
ambient plasma at an initial temperature of 1 GeV or 10 GeV respectively. 

Since our formalism above is not very rigorous we shall also estimate below the asym- 
metry in the diffusion picture ignoring any plasma motion due to pressure gradients. For 
a black hole of initial mass 4 x 10 11 Mp; the volume fraction Vribh, where V = |vrrf 00 
is 2 x 10 -11 . Taking into account the entropy factor, the asymmetry is 2 x 10~ 5 of the 
asymmetry in standard electroweak baryogenesis models. For a black hole of initial mass 
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8 x 10 l0 Mpi the volume fraction is 4 x 10~ 7 . So the asymmetry is about 4 x 10 -4 of the 
asymmetry in standard electroweak baryogenesis models. 

In the diffusion picture we have ignored the fact that once the black hole has evaporated 
completely energy from hotter regions will diffuse out and heat up cooler regions. This 
will effectively increase rioo- Furthermore, if we relax the constraint that the black holes 
do not dominate the Universe until they decay this may enhance the final asymmetry. 

We point out an interesting phenomenon that can occur in the regions surrounding 
the black hole where the temperature rises but does not get much higher than T ew . In 
these regions the effective potential changes from the zero temperature potential to a high 
temperature potential with a local minimum at < >= v as shown in fig. 1. In these 
regions the Universe is trapped in a state with < >= v and tunnels to the < >= 
vacuum by bubble nucleation. Thus in these regions bubbles of < >= are created. As 
long as the < <p >= v vacuum is higher than the < >= vacuum these bubbles expand, 
as in fig. 2a. However as this region cools the two vacua become degenerate and then the 
state with < >= v becomes the true vacuum. If the < >= bubbles have not already 
disappeared after collisions (leaving the Universe in the < >= vacuum state) then 
the < >= bubbles now begin to shrink, as in fig. 2b. Baryons are now created in the 
inner regions of the shrinking walls. Furthermore, < <p >= v bubbles may nucleate inside 
the < >= bubbles (depending on the critical bubble size). As these expand baryons 
are created in their walls. The above phenomenon may occur in the diffusion picture for 
r > r 10 o in the transition phase after the black hole evaporates away and energy from 
hotter regions moves out and heats up colder regions. In the case of heat transfer through 
bulk motion of the plasma this mechanism may occur at some distances away from the 
black hole in the early stages before the steady state is achieved. 

We now consider a second order phase transition. In the hot regions, the symmetry 
is restored and the field settles at the minimum of the potential at = 0. As the region 
cools below T ew the minimum shifts from = 0. < >, in this case, is not due to Hubble 
expansion but is due to the cooling of the plasma as the plasma expands. Hence, unlike in 
the standard cosmology, < > is large. The non-zero < > acts as a chemical potential 
for baryon number during the phase transition and a net baryon number is created by 
B+L violating sphaleron processes. The baryon asymmetry generated in a second order 
electroweak phase transition is given by H3[] 

n B /s « 10- 19 e[-T/T] 100 /# 100 (28) 

where H wo is the value of the Hubble constant when the temperature of the Universe 
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is 100 GeV, T/T is in the frame of the plasma when the temperature has fallen to 100 
GeVand e is a measure of CP violation. 

To estimate T/T in the rest frame of the plasma when the temperature has fallen to 
100 GeVwe note that at r 10 o in time dt a shell of thickness dr = Vio^dt moves out at a 



temperature Tioo- As the shell moves out T = f^r^r;^ ~ —Tvy/r, since eqs. fl22"D and 



( 2lD imply that T is approximately proportional to r . dt' is the time differential in the 



plasma frame. Therefore [— T/T] wo = t>ioo7ioo/ r ioo- 

Using values of ri o and t>ioo obtained earlier we get [— T/T] wo = 1 x 10 _1 7^ 00 GeV and 
3 x 10 _2 7 2 00 GeV for Tu equal to 1 GeVand 10 GeVrespectively. In the standard scenario 
[— T/T] ioo = -^ioo = 10~ 14 GeV. Therefore the baryon asymmetry that we obtain in this 
scenario, taking into account the volume fraction and the entropy factor, is 

n B /s = 3 x 10~ 9 7iooe GeV" 1 (29) 

for Tjj = 1 GeV and is 

n B /s = 9 x 10" 9 7i 00 e GeV" 1 (30) 

for Tu = 10 GeV. Since 7100 ~ 10 3 for our choices of rj clearly one can obtain a sufficiently 
large baryon asymmetry in our scenario with a second order phase transition. It is inter- 
esting that while the asymmetry generated in a first order electroweak phase transition 
is suppressed by the Lorentz factor, the asymmetry is enhanced by the Lorentz factor for 
a second order phase transition. Also note that even if viscous effects reduce the value 
°f 7ioo to 1 we can obtain sufficient asymmetry for e ~ 10" 1 — 10~ 2 (for t> 10 o not much 
smaller than 1). Further note that the possibility of creating a sufficient asymmetry is not 
very sensitive to the value of 77 as long as T c is large enough (^ 200 GeV) to ensure that 
fioo is not much smaller than 1. 

Since the distance between black holes is much greater than rjj, the distance at which 
the outward moving plasma attains the ambient temperature, the heating of the plasma 
by one black hole does not directly affect the plasma around neighboring black holes. 
Once the outward moving plasma reaches Tu it merges with the ambient plasma. As the 
energy density of the black holes does not dominate the energy density of the universe the 
heating of the ambient plasma between black holes, which could affect the temperature 
profile of neighboring black holes, is also negligible. Finally, baryons carried out by the 
outward moving plasma into the ambient plasma can subsequently enter the hot region 
around another black hole due to random Brownian motion across the intermediate region. 
However as the baryons cross rjj of another black hole the plasma moving out from the 
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second black hole would push them out thereby preventing them from reaching rioo of the 
second black hole, where sphaleron processes in equilibrium could destroy them. 

We now briefly comment on the possibility of creating hot regions by the decay of 
massive particles. Massive particle decays will lead to the development of a collimated jet 
of particles, with each particle in the jet contributing to the heating of the plasma. From 
eqn. (7) for the energy loss of a particle created by the decay of this massive particle it seems 
very difficult to get sufficient energy to heat the region to electroweak scale. However, note 
that the expression for the energy loss as given in eqn. (7) has been derived by neglecting 
very hard collisions. Inclusion of hard collisions may lead to a rapid branching of the 
initial decay product which may significantly increase the energy loss. Also note that 
baryogenesis will only take place in regions that are at least as wide as the sphaleron size. 



IV. CONSTRAINTS ON THE HIGGS MASS 



In the electroweak baryogenesis scenario in the context of the Standard Model if the 
high temperature vev in the broken phase is small then the sphaleron processes are still 
active. These sphaleron processes will destroy the baryon asymmetry created during the 
phase transition. A large enough high temperature Higgs vev in the broken phase to avoid 
the wash-out of the asymmetry implies a small Higgs mass which is tightly constrained by 
the experimental lower bound on the Higgs mass. In our scenario, on the other hand, the 
requirement of a small Higgs mass does not apply, as we show below. Of course, this is 
relevant only if viscous forces are sufficient to decrease the suppression of the asymmetry 
due to the Lorentz factor for a first order electroweak phase transition. A large Higgs 
mass also implies that the electroweak phase transition is likely to be of second order 
in which case the asymmetry generated at the electroweak phase transition in standard 
electroweak baryogenesis is very small ||43|| . As we have already indicated, we generate 
sufficient asymmetry with a second order phase transition in our scenario. 

The baryon number density after a first order electroweak phase transition is given by 

n B = n Bi exp(— / Tdt) = n Bi exp(— / T/TdT) (31) 

J J100 

where n B i is the baryon number density created during the phase transition, T(T) is the 
rate per unit time of B violation and we choose the upper limit T' to be Tu, the ambient 
temperature of the universe. More precisely, one could put T' as the temperature at which 
the B violation rate becomes less than the Hubble expansion rate but the above choice is 
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sufficient. For T we use —(jv/r)T and r = r 100 T 100 /T. 7 is given by eqn.(^TJ) and we let 
v = l. 

The perturbative rate of baryon number violation in the broken phase is given by |44 



T = yfotwT 



M 7 W 
a 7 w T 7 



exp[-E sph /T] (32) 



where E sph = 2M w (T)/a w , with M W (T) = M W (0)^J1 - T 2 /T* w and T ew ~ 100 GeV is 
the critical temperature. The above is valid for 2M W (T) C T < 2M w (T)/aw, i-e., it 
breaks down close to T ew . y is obtained by comparing the B violation rate to the sphaleron 
rate per unit volume in eqns. (2.17) and (2.14) of ref. ||44|| . Choosing B-L=0 and three 



families we get y=20. / includes other factors found in eq. (2.19) of ref. [44 



The non-perturbative lattice calculation of the sphaleron rate per unit volume in the 
broken phase close to the critical temperature T ew is given by fl45" 

Fsph/V = xa^T 4 , (33) 

where x varies from 10~ 8 to 10™ 5 for different values of the couplings in the theory. 
Therefore T(T) is given by 

T = xyafyT . (34) 

We can model the rate of baryon number violation for the range T = 100 GeV — TJy as 

T(T) = xya 4 w Texp[-E sph /T}T 7 ew /T 7 (35) 

where we ignore the temperature dependence in M w . This should hence give us an upper 
bound on the actual rate. 

We then find that the reduction in the baryon number density due to sphaleron pro- 
cesses is negligible. (The fraction of baryons destroyed before the temperature of the 
expanding plasma becomes Tu is 10~ 10 and 10 -11 for T v equal to 1 GeVand 10 GeV. re- 
spectively.) This is true even for much larger values of x, including the larger perturbative 
estimates of x shown in fig. 4 of ref. |[45|| . 

The expression for the asymmetry in the second order phase transition case, adapted 
from the results of ref. jKJ, takes into account any destruction of the asymmetry when 
the value of the asymmetry at any time becomes greater than the equilibrium value. 

The above arguments, coupled with the observation that our model works for a second 
order electroweak phase transition, implies that the constraint of a lower bound on the 
Higgs mass does not apply to our scenario. However, as we have pointed out earlier, 
insufficient CP violation in the Standard Model requires one to consider extensions of the 
Standard Model. 
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V. CONCLUSION 



The baryon asymmetry that is created in our scenario is inhomogeneously distributed 
in the universe. However an estimate of the distance traversed by the baryons beyond rjj 
(in the plasma motion case) due to random Brownian motion indicates that the baryons 
produced in our scenario homogenize themselves on time scales of the order of the lifetime 
of the black holes. Particles traveling through a plasma suffer collisions and the actual 
distance traveled by a particle in the time r is \f~Nl p , where N is the number of collisions 
and l p is the mean free path |L6f. Now N = r/(l p /c). Therefore y/Nl p = (tcZ p )s. The 
mean free path=l/(em) and, as stated earlier, an ~ 0.02T. If we substitute the black 
hole lifetime for r then the distance traversed by the baryons beyond rjj due to random 
Brownian motion is 5 x 10 9 GeV -1 and 1 x 10 s GeV -1 for ambient temperatures of 1 
GeVand 10 GeVrespectively. Since the distance between black holes is 5 x 10 9 GeV -1 and 
1 x 10 8 GeV -1 for ambient temperatures of 1 GeVand 10 GeVrespectively and the inverse 
Hubble parameter is of the order of the black hole lifetime this indicates that the baryon 
density inhomogeneities can be wiped out on time scales of the order of few times the 
black hole lifetime. 

Neutrino inflation and baryon diffusion can also homogenize any baryon number den- 
sity fluctuation prior to the onset of primordial nucleosynthesis. Neutrino inflation is the 
process of heating up of overdense baryon density lumps in the early Universe by neutrinos 
which causes the lumps to expand, or inflate, thereby reducing their overdensity. Besides 
this mechanism, neutrons and protons also diffuse out of overdense baryon density lumps. 
Mechanisms of diluting overdense baryon density lumps are discussed in ref. ||47|| . The 
above arguments also apply to the diffusion scenario where the size of the baryon density 
inhomogeneity is given by r WQ . 

As we mentioned earlier, the scenario in ref. [18] has similar features to our work. 
In this paper one obtains the temperature profile outside the black hole in a manner 
similar to that used for stellar interiors. As we have pointed out, the assumption that 
heat transfer is primarily due to diffusion, as in stellar interiors, is not necessarily valid 
for the plasma surrounding a primordial black hole. 

There are many open questions which remain to be explored. The issue of black hole 
formation with the required number density and mass spectrum for our scenario requires 
to be discussed in the context of some realistic model such as a first order phase transition 
occurring earlier in the Universe. Also one must keep in mind observational constraints 
on massive black holes [|8-p| which may lead to strong restrictions on the spectrum of 
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masses of primordial black holes. A more complete numerical calculation can give us an 
estimate of the temperature of the plasma close to the black hole and fix 77. Finally viscous 
effects can be included in the relativistic Euler equation to more accurately ascertain the 
velocity of the plasma. 
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FIGURE CAPTIONS 

Fig.l: Effective potential when the region superheats with < <p >= v becoming 
metastable. 

Fig.2: (a) Nucleation of < >= bubble in the metastable phase with < >= v. 
(b) Nucleation of < <fi >= v bubble inside a shrinking < <fi >= bubble. 
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